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The ground state energy and the lowest excitations of a two dimensional Wigner crystal in 
a perpendicular magnetic field with one and two electrons per cell is investigated. In case of two 
electrons per lattice site, the interaction of the electrons within each cell is taken into account exactly 
(including exchange and correlation effects) , and the interaction between the cells is in second order 
(dipole) van der Waals approximation. No further approximations are made, in particular Landau 
level mixing and incomplete spin polarization are accounted for. Therefore, our calculation comprises 
a, roughly speaking, complementary description of the bubble phase (in the special case of one and 
two electrons per bubble), which was proposed by Koulakov, Fogler and Shklovskii on the basis of a 
Hartree Fock calculation. The phase diagram shows that in GaAs the paired phase is energetically 
more favorable than the single electron phase for, roughly speaking, filling factor / larger than 0.3 
and density parameter Vs smaller than 19 effective Bohr radii (for a more precise statement see Fig.s 
4 and 5). If we start within the paired phase and increase magnetic field or decrease density, the 
pairs first undergo some singlet- triplet transitions before they break. 

PACS: 71.10.Hf Non Fermi liquid ground states, 

71. 10. Li Excited states and pairing interactions in model systems, 

73.40.Hm Fractional Quantum Hall Effect 

I. INTRODUCTION 

The Wigner crystal (WC) and the Fermi liquid are states of matter which are established already for decades. Later 
on, the Laughlin liquid as a special state in high magnetic fields has been added. In both liquid states, wave-functions 
of the electrons have an appreciable overlap. In the extreme Wigner crystal limit this overlap can be neglected. The 
issue of this paper is if there is a state between crystal and liquid, where a finite number of particles overlap as in 
the liquid, but interacts with the neighboring cluster as in the Wigner crystal (paired or clustered Wigner crystal). 
The charge density-wave (CDW) state, mostly investigated in Hartree Fock (HF) approximation, is a alternative 
intermediate phascErQ In the latter case, there is a density modulation, but the density is so high that overlap plays 
an important role. (Observe that in some papers the CDW state is also called WC.) What we call here 'paired Wigner 
crystal' can be viewed as a derivative of three other states. It is as (i) a Wigner crystaH with two electroirs-per cell, 
(ii) a crystallized paired electron liquid known from superconductivity and theories for the / = § stateETu , (iii) a 
state of liquid islands, if the clusters are large. 

The first paper which considered the possibility of a paired WC in 3 dimensions (3D) without a magnetic field is 
Rcf. ^ using a cell approximation for the inter-cell interaction and a variational ansatz for the intra-cell problem. In 
a previous paperU it was shown, that the model of Rcf. |l^ can be solved exactly, even if fluctuation corrections in 
second order van der Waals approximation are included additionally. It turned out that in 3D without a magnetic 
field, the paired WC is energetically lower than the conventional WC for Vg < 9, but in this density region the 
Fermi liquid is probably already lower than-bpth WCs. More recently, this issue has been investigated in 2D with a 
magnetic field using the HF approximation.ll3 The result is that in higher Landau levels (LL) the ground states_ai£ 
phases with circular or stripy clusters. Meanwhile, numerical solutions for finite size systems have been publishectjO, 
which confirm these findings qualitatively within their model. All previous work on this issue neglect LL mixing and 
incomplete spin polarization. In the present approach, we investigated clusters with one and two electrons. In the 
latter case, the interaction within each cluster is exactly accounted for (including exchange and correlations). In both 
cases, the inter-cell interaction is considered in second order van der Waals approximation (no exchange, Coulomb 
interaction in dipole approximation). In Sect. II we describe the solution of the Schrodinger equation in more detail 
and in Sect. Ill our results are presented. Sect. IV gives a summary and an overview. 
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II. SOLUTION OF THE SCHRODINGER EQUATION 



The Hamiltonian of a system of N electrons at each lattice site and a positive compensating background of density 
no reads (atomic units h = m = e = 1 are adopted) 



with the kinetic energy 



2m* 



n k 



Pnfe + -A(u„fe) 
C 



the intra-cell e-e-interaction 



the inter-cell e-e-interaction 



VI 
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n k^k' 



the electron- background interaction 
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In these definitions, m* is the effective mass, /3 the inverse background dielectric constant, R° are the lattice sites, and 
r„fe = R° + Unk the electron coordinates. We use the symmetric gauge A = ^ B x u. Latin letters n and n' denote 
the lattice sites and k and k' the electrons within a cell (or cluster). If the distance between the cell centers is large 
compared with the average distance between the electrons within a cell, the inter-cell interaction can be expanded in 
a multi-pole series up to second order 



n k 



and the background- background interaction 
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with the zero and second order terms 
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the force constant matrix 



the dipole tensor 
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Cnk,n'k' = +pN^To for ink) = {n'k') 
= -pNTnn' else 
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and the definitions To = X^„^o "^O" ~ ^n ~ ^n' ■ rearrangement of terms we obtain the following form 



2 



H = Efj^^ + (T + T/-,, J + V::!^ (13) 
with the Madelung energy for a lattice with N electrons at each lattice site 

By introducing the center of mass (cm.) coordinate U„ and some relative coordinates within each cell using a 
orthogonal coordinate transformation, the last 3 terms of (|l^ ) can be decomposed into a cm. Hamiltonian Hem. and 
a lattice sum of internal (relative coordinate) Hamiltonians Hrei = Tim Hrei,n, where Hrei,n contains only the {N — 1) 
relative coordinates in cell n. Thus, only Hem. shows a coupling between the cells. This fact is completely analogous 
to the treatment of interacting quantum dot lattices considered in Ref. ^ Consequently, the total Hamiltonian ( p^ ) 
decomposes into 3 contributions 

H = EfJ""'^ + Hem + Hrel (15) 

which are considered in turn. 



A. Madelung energy 



Generally, we obtain from the conventional Madelung energy {N = 1) with the same lattice constant by 

jT^Mad _ j\j2jTjMad^ Here, howcvcr, we compare phases with the same mean density. In the hexagonal lattice, the 



Madelung energy per electron ey} for a given density parameter rs = l/y^Tmo can be deduced from the data in Ref. 
I providing e^J"'^ = -1.106103 /3 N^^^ rj^. 



B. Center— of-mass energy 

The cm. part of the Hamiltonian reads 



N 

P„ + -A(U„) 

c 
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with 



u .c ,.u 



C„,„' = +/37V2To for n^n' 
= -/?7VT„„, else 



After the usual phonon transformation 
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where N^. is the number of cells, the Hamiltonian decouples with respect to the lattice sums and we obtain 

BZ 



Hem — ^ ^ 



(21) 



with the magneto-phonon Hamiltonian 
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N 



N 



A(u;) 



and the dynamical matrix 



(22) 



(23) 



The eigenvalues of Hem and //q depend on N only through the dependence of C. The explicit N dependence in 
( p^ ) and ( p2| ) cancels. The last statement is obvious if we consider that the explicit N can be removed by rescaling 
the displacement vector U by a factor of ^/N. 
The eigenvalues of (122) are e.g. given in Ref. 



E{n+, n_) = (n+ + i) w+ + (n_ + i) ; n± = 0, 1, 2, . 



This provides a ground state energy 



BZ 



(24) 



(25) 



where 



UJ± 




(26) 



with uJq 



^(Cii + C22), A = Cii — C22, = is the cyclotron frequency with the effective mass, and Ci_k are 
the Cartesian components of Cq. If we define an auxiliary tensor S through Cq — p Sq with p — 2l3N/a^ (a is the 
lattice constant), then Sq depends only from the geometry of the lattice, i.e. it is a fixed tensor for the hexagonal and 
cubic lattice and depends only from the b/a ratio for a rectangular lattice. The explicit form of Sq is complicated 
and not of common interest. The sums involved have to be done numerically anyway. The interaction parameter p 
for the hexagonal lattice reads in terms of the more convenient parameters rs and uq 




3/2 



3/2/V-1/2 



(27) 



Because of the optical selection rules, ijJ± for q = agree with the excitation energies AS for long wavelength 
(infrared) radiation. 

It can be easily seen that for B = 0, ijJ± ~ eigenvalues{Cq) — y/p ^ eigenvalues{Sq) ■ Because of the universal 
character of Sq, the excitation energies in units of ^/ p comprise all data on the energies for hexagonal lattices in 
this limit (see Fig.l). In connection with formula ( p7[ ) this means that the phonon frequencies of the paired WC, as 
compared with the single electron case of the same electron density, are smaller by a factor of 1/^/2. Additionally, 
the frequencies (and the cm. contribution to the total energy) for the same N decay with the density parameter like 
-3/2 
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For finite B, an universal function comprising all data cannot be defined. However, if we use co* as energy unit and 
uj*'^ as unit for the interaction parameter p, then the result depends only on p and shows no explicit dependence on 
B. 

For large i?, ( p6| ) gives 0;+ — lo* + trace{Cq) /2lu* and aj_ = determinant{Cq) /lo* (in agreement with Ref. and in 
the units mentioned above, this reads uj^[uj*.] =^ 1 • trace{Sq)/2 and — p[{uj*)'^] ■ deter'minant(Sq) . 

Thus, the results are determined by the interaction parameter and two universal functions, namely trace{S) and 
determinantlS) , which are shown in Fig. 2. In the limit B 00 we obtain the result of noninteracting electrons 
— LO* and lo^ =0, as to be expected. 




(qx.qy) 

FIG. 2. The functions trace(Sq) and deter minant{Sq) determining the excitation frequencies in the high field limit. 
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Direct comparison of excitation energies with experiments is not the aim of this work. What matters here is only 
the contribution of the 'zero point vibrations' to the total energy. 



C. Internal energy 



For two electrons per cell, the internal (relative coordinate) energy can be calculated easily. (Larger N call for a 
larger numerical effort and are postponed to a succeeding work.) If we introduce apart from the cm. the relative 
coordinate u„ 



U„ = U„2 — U„i 

then all relative Hamiltonians read (the index 'n' is omitted) 



'.I 







P + ^A(u) 


[2m* 



1 X. P 

- u D u+ — 

2 2u 



(28) 
(29) 



(30) 



where p = — iVu and D = ^Tq = jpSgl with 5*0 being a lattice sum of the same type as the components of Sik- (Do 
not mix up the interaction parameter p and the vector of the momentum operator p.) It should be emphasized that 
the electronic mtercell interaction contributes a harmonic term to the miracell problem to Hrei, namely the last but 
one in (refH-rel). As to the calculation of the eigenvalues of Hamiltonians of the form (^ ) we refer to Ref. |l^. 
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III. RESULTS 



If not otherwise indicated, all results and parameters are given in effective atomic units for GaAs {m* = 0.067, f3 = 
1/12), .i.e., energies (and frequencies) in la.u.* = 4.65 • 10^'* double Rydberg = 12.64 meV and length in la.u* ~ 
1.791 • lO^Bohr = 0.9477 • lO^A. (If the units of a quantity is indicated explicitly, then it is given in brackets after 
the quantity.) Observe that the effective density parameter in current typical experiments is of order 2. 



A. Contributions to the total energy 



For larger r^, the energy difference between both phases is tiny (see Fig. 3). As an example, at r, = 10 (at the 
right boundary of Fig. 3) it amounts to 2 • 10^^, calling for high precision computations (of order lO^'^a.ii.*). This 
fact is amazing because the Madelung energy and the cm. energy are both very different (see Fig. 3). This difference 
is compensated almost completely by the relative energy of the paired phase. 

1 I ^ ^ ^ ^ , ^ ^ ^ ^ 1 




FIG. 3. Comparison of the contributions to the total energy in the single electron (N=l) and the paired (N=2) Wigner 
crystal for zero magnetic field. 



B. Phase diagram for single electron versus paired Wigner crystals 

Fig.4 shows the phase diagram for the single electron versus the paired WC phase in the rg- B plane. In shaded 
regions, the electron pairs arc in singlet state and in triplet states elsewhere. The spin configuration of the paired 
phase has been determined irrespective of the question, if the single electron phase is below both paired (singlet or 
triplet) phases. The thin full lines separate regions with different relative (internal) angular momentum m of the 
pairs, which has the lowest energy. The value of m is indicated in each region as well. (Even m belongs to the state 
and odd m to the triplet state.) The thick full line separates the regions where the single electron WC is the ground 
state from regions, where the paired WC has the lower energy. We want to emphasize that this figure does not say 
anything about the issue if any other (e.g. liquid) phase is lower in energy than the phases considered here. The kinks 
in the phase boundary are connected with changes of the total spin of the paired phase. 
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rs [a-u. ] 




Ts [a-u- ] 

FIG. 4. Phase diagram in the r^- B plane (see text for details) a) for small B and b) for small r, 



First, we discuss the case B = 0. For above 19, the energetically lower state is the single electron WC and below 
it is the paired phase. Clearly, in the limit of high densities (small Vs) any liquid state should win. It is physically 
imaginable, that in increasing the density, the state starts with the single electron WC, then adopts the paired state 
before it goes over to the liquid. If there are states with 3, 4, ... electrons involved, will be the focus of further 
investigations. 

For a certain finite B, the density, where the transition to the paired ground state occurs, becomes higher. In other 
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words, for higher B one needs less dilute electron systems to produce a single electron WC. This is physically obvious, 
because the magnetic field helps to localize the electrons. In other words, lowering the density cases pair breaking. 
Now we discuss the pair breaking process as a function of B for fixed r^, say rg = 2. If we start with B = and 
increase B, the paired state undergoes first some spin transitions from singlet to triplet and vice versa before the pairs 
break in the triplet state. It should be mentioned that pair breaking and spin transition are not coupled together. 
The reason is that the Zeeman energy is too small to play a decisive role. 



N-P 

1 N — ^ 

m=0 


/ 


m=-1 


N=1 


m^^3z_l!!lZ^ — 



10 . 20 



r3 [a.u.] 

FIG. 5. Phase diagram in the r,- f plane (see text for details). 



Fig. 5 shows the the phase diagram again, but the magnetic field axis is replaced by the filling factor / — 
14.64/(i3[T] r^[a.u.*]) (using the parameters of GaAs). This plot has the advantage of showing the physically impor- 
tant parameter / directly, but does not allow the discussion of the limit B = 0, because it is transformed to / — > oo.) 
Again, thin lines separate regions with different m (and therefore total spin) of the paired phase and the thick line 
separates the single electron WC ground state from the paired WC. As to be expected, the single electron WC is 
restricted to the low density- small filling factor region. It is clear that in a complete phase diagram, the region, 
which is here attributed to the paired phase, can be further restricted by the inclusion of all kinds of liquid phases 
(and possibly WC phases with N=2,3,... electrons per cell). 

C. Comparison of Wigner crystal phases with Laughlin liquid 

Because the Laughlin liquid is restricted to discrete filling factors, we compare in Fig. 6 the energy (or, more 
conveniently, ^^(e — lUc/2)) of different phases for fixed / = 1/3 as a function of rj. The Zeeman energy is omitted in 
this figure because the paired phase is triplet over the whole Ts range. Therefore it would have contributed for a given 

the same amount to all phases, but introduce an unpleasant pole at = into the curves. (In all calculations, 
however, the Zeeman energy has to be taken into account because of the different spin configurations in the paired 
WC.) The curve for the Laughlin liquid from Ref. |l^ includes Landau level mixing (as our results do) by adding a 
variational correlation factor to the Laughlin function. (The value at = corresponds to the result of the original 
Laughlin function.) As seen in Fig. 6, the Laughlin liquid and the single electron WC are the ground state for rs 
below and above 15, respectively. This is in qualitative agreement with experiment. In our result for / = 1/5 (not 
shown), the single electron WC is the ground state over the whole range of densities. This is not quite correct, because 
experiments suggest that the domain of the WC begins at / = 1/7. Inclusion of a variatiapal correlation factor to the 
WC wave-function fixes the phase boundary between the Laughlin liquid and the WCcjEj. This amendment could 
not be added here, because it would have become too complicated for the paired WC and in comparing both WC 
phases we needed a description on the same footing. 
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IV. SUMMARY AND OVERVIEW 



We have shown that the paired WC is energetically more favorable than the single electron WC for higher densities 
and lower magnetic fields. If we start with the paired WC, lowering the density and increasing the magnetic field 
will break the pairs. At least in GaAs, the Zeeman energy plays a less important role than the angular momentum 
dependence of the internal (relative coordinate) energy. Observe that changing the internal angular momentum (from 
even to odd) is coupled with a change of the total spin of each cluster due to the Pauli principle. Therefore, the 
paired phase undergoes first some singlet- triplet transitions before the pairs break and the single electron WC is 
reached. The singlet- triplet transitions have the same physical origin as the corresponding transitions in quantum 
dots, because they are driven by the internal energy, which is similar in either system. 

In our approach, the interaction within each cluster is exactly accounted for (including exchange and correlations). 
In both cases, the interceW interaction is considered in second order van der Waals approximation (no exchange. 
Coulomb interaction in dipole approximation). The validity of the second order approximation has been checked in 
Ref. |l^ by estimating the third order contribution, which amounts to a few %. We did not include higher order 
contributions from two reasons. First, they spoil the exact decoupling and make the paired WC intractable. Second, 
these corrections have most likely the same sign in the single electron and the paired phase and they cancel partly in 
comparison of both energies. The neglect of inter-cell exchange (but keeping the Coulomb correlations) is justified in 
the low density and high magnetic field limit, because the exchange energy drops off exponentially, but the correlation 



energy like l/r. It seems to be worthwhile to elaborate on the differences of this work and Ref. 12. They targeted 
on high LLs (low magnetic fields), and high electron densities {N^l >> 'r~^ >> 1). LL mixing and incomplete spin 
polarization are ruled out and the completely filled LLs are frozen out and influence the upper LL only through a 
screening in the e- e- interaction. Therefore we called both treatments in the abstract 'complementary'. What seems 
to be interesting is that despite approaching reality with different assumptions and virtually from opposite sides, the 
effect of clusterization of the liquids is found in both treatments. 
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